UNCLASSIFIED 


FOR  REPRODUCTION  PURPOSES 


AD-A211  180 


PORT  DOCUMENTATION  PAGE 


CTE 


2b.  DECLASSIFICATION  /  DOWNGRA 


4.  PERFORMING  ORGANIZATION  R 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

JAI  Associates,  Inc. 

6b.  OFFICE  SYMBOL 
(If  applicable) 

6c  AODRESS  (City.  State,  and  ZIP  Code) 

P.  0.  Box  293,  Mountain 

View,  GA  94  OL 

8a.  NAME  OF  FUNDING  /SPONSORING 
ORGANIZATION 

U.  S.  Army  Research  Office 

8b.  OFFICE  SYMBOL 

Of  applicable) 

8c  ADDRESS  (City.  State,  and  ZIP  Code) 

P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 

|  3  .  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

0  Si  S~b"7S~  $'€6-  S _ 


7a.  NAME  OF  MONITORING  ORGANIZATION 
U.  S.  Army  Research  Office 


7b.  ADDRESS  (City,  Sfata,  and  ZIP  Code) 

2  P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  I  PROJECT 
ELEMENT  NO.  I  NO 


11.  TITLE  (Include  Security  Classification) 

Unsteady  Interaction  of  a  Rotor  with  a  Vortex,  AIAA  Paper  AIAA-89-1848 


12  PERSONAL  AUTHOR^) 

G.  R.  Srinivasan  £  W.  J.  McCroskey 


13b.  TIME  COVERED 
FROM _ TO 


14.  DATE  OF  REPORT  (Tear.  Month,  Day)  IlS.  PAGE  COUNT 


13a.  TYPE  OF  REPORT 

Reprint 


16.  SUPPLEMENTARY  NOTATION  ,  ,  .  ,  .  , , 

The  view,  opinions  and/or  findings  contained  in  this  report  are  those 

of  the  authqr(s) .and  should  not  be  construed  as  an  official  Department  of  the  Army  position, 

nnlirv  r\r-  Hon  cion  nnloco  on  Hbci  onafod  Ktr  nt-nor  HnnnmonP  of  ihn  J 


18.  SUBJECT  TERMS  ( Continue  on  reverse  if  necesssry  and  identify  by  block  number) 

Unsteady;  Euler  £  Navier*-Stokes ,  Parallel  and 

Oblique  Blade-Vortex  Interaction! 


COSATI  COOES 


GROUP  SUB-GROUP 


’9  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 


The  unsteady,  three-dimensional  flow  field  of  a  helicopter  rotor  blade  encountering  a  passing 
vortex  is  calculated  by  solving  the  Euler/thin  layer  Navier-Stokes  equations  by  a  finite-difference 
numerical  procedure.  A  prescribed  vortex  method  is  adopted  to  preserve  the  structure  of  the 
interacting  vortex.  The  cases  considered  for  computation  correspond  to  the  experimental  model 
rotor  test  conditions  erf  Caradonna  et  al.  and  consist  of  parallel  and  oblique  interactions.  Com¬ 
parison  of  the  numerical  results  with  test  data  show  good  agreement  for  both  parallel  and  oblique 
interactions  at  subsonic  and  transonic  tif 


89 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 

□  UNCLASSIFIEDAJNLIMITEO  □  SAME  AS  RPT  □  OTICUSERSl  Unclassirieu 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL  |22b.  TELEPHONE  (Include  Area  Code)  22c.  OFFICE  SYMBOL 


DD  FORM  1 473, 84  mar 


S3  APR  edition  may  b«  used  until  exhausted. 
All  other  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 
UNCLASSIFIED 


Al  AA-89-1 848 

UNSTEADY  INTERACTION 
OF  A  ROTOR  WITH  A  VORTEX 

G.  R.  Srinivasan, 

JAI  Associates,  Inc.,  Sunnyvale,  CA 
W.  J.  McCroskey, 

U.  S.  Army  Aeroflightdynamics  Directorate-AVSCOM 
NASA  Ames  Research  Center,  Moffett  Field,  CA 

Accesior.  For 

jNTJS  CRAW 
OTIC  TAB 
Unannounced 
Justification 

By _ _ 

Oist- (button/ 


'  Avail 
Dist  ;  so 

/Mao 

AIAA  20th  Fluid  Dynamics,  Plasma  Dynamics 

and  Lasers  Conference 

Buffalo,  New  York  /  June  12-14,  1989 


UNSTEADY  INTERACTION  OF  A  ROTOR 
WITH  A  VORTEX 

G.  R.  Srinivasan  * 

JAI  Associates  Inc.,  Kifer  Court,  Sunnyvale,  California 
W.  J.  McCroskey  ** 

U.  S.  Army  Aeroflightdynamics  Directorate  -AVSCOM 
NASA  Ames  Research  Center,  Moffett  Field,  California 


Abstract 

The  unsteady,  three-dimensional  flow  field  of  a  heli¬ 
copter  rotor  blade  encountering  a  passirg  vuriex  is  calcu¬ 
lated  by  solving  the  Euler/thin  layer  Navier-Stokes  equa¬ 
tions  by  a  finite-difference  numerical  procedure.  A  pre¬ 
scribed  vortex  method  is  adopted  to  preserve  the  structure 
of  the  interacting  vortex.  The  cases  considered  for  com¬ 
putation  correspond  to  the  experimental  model  rotor  test 
conditions  of  Caradonna  et  al.J,J  and  consist  of  parallel  and 
oblique  interactions.  Comparison  of  the  numerical  results 
with  test  data  show  good  agreement  for  both  parallel  and 
oblique  interactions  at  subsonic  and  transonic  tip  speeds. 

Nomenclature 

a  =  speed  of  sound 

a0  =  vortex  core  radius,  see  Eq.  (6) 

doo  =  free  stream  sound  speed 

C  =  characteristic  length  scale,  chord  of  the  blade 

Cw  =  chord  of  the  vortex  generating  wing 
Cp  =  pressure  coefficient,  based  on  local  dynamic 

pressure 

e  =  total  energy  per  unit  volume 

F,G,H  =  flux  vectors 

J  =  Jacobian  of  the  coordinate  transformation 

Moo  =  free  stream  Mach  number 

Mtip  =  tip  Mach  number  of  the  rotor  blade 

p  =  static  pressure 

Poo  =  free  stream  pressure 

Pr  =  Prandtl  number 

Q  =  flowfield  vector 

r  =  radial  distance  from  the  vortex  center 

rg  =  rotor  reference  station  normalized  by  R 

R  =  rotor  radius 

R(i)  —  rotational  matrix,  see  Eq.  (4) 

Re  -  Reynolds  n>’rr>b'- 

S  —  viscous  flux  vector 

«oo  =  free  stream  velocity 

u,v,w  =  velocity  components 

*  Senior  Research  Scientist,  Associate  Fellow,  AIAA. 

**  Senior  Staff  Scientist,  Associate  Fellow,  AIAA. 


ve  —  tangential  velocity  of  vortex 

U,  V,  W  =  contravariant  velocity  components 
x,y,z,t  =  inertial  coordinates 

xa,  z0  =  vortex  offset  position  relative  to  rotor  axis 
xv,zv  =  distance  from  blade  leading  edge  to  vortex 
z,y,z,t  =  blade  attached  coordinates 
7  =  ratio  of  specific  heats 

T  =  circulation  for  the  vortex 

T  =  dimensionless  strength  of  the  vortex,  r/(aODC) 

p  =  viscosity  coefficient,  advance  ratio 

£  i  Vi  C> T  =  generalized  curvilinear  coordinates 

Poo  -  free  stream  density 

p  =  density 

i/>  =  azimuth  angle 

Cl  =  angular  velocity  of  the  rotor 

Introduction 

The  accurate  simulation  of  the  flowfield  of  a  heli¬ 
copter  rotor  is  still  one  of  the  most  complex  and  challenging 
problems  of  applied  aerodynamics.  This  is  true  in  spite  of 
the  availability  of  the  present  day  supercomputers  of  Cray- 
2  class  and  improved  numerical  algorithms.  The  flowfield 
of  a  rotor  in  forward  flight  is  highly  three-dimensional,  un¬ 
steady  and  viscous,  with  pockets  of  transonic  flow  near  the 
blade  tips  on  the  advancing  side  and  regions  of  dynamic 
stall  on  the  retreating  blades.  In  addition,  the  blades  also 
shed  complex  vortical  wakes;  the  concentrated  tip  vortex 
of  each  blade  generally  passes  close  to  the  following  blades. 
The  close  encounter  of  these  (force  free)  concentrated  vor¬ 
tices  with  the  rotor  blades  is  often  the  cause  of  unsteady 
load  fluctuations  and  acoustic  noise.  Some  of  these  com¬ 
plexities  associated  with  the  flowfield  are  delineated  in  a 
schematic  of  helicopter  rotor  in  forward  flight  shown  in 
Fig.  1. 

The  spiralling  vortex  sheet  emanating  from  each  of 
the  blades  of  the  rotor  has  a  profound  influence  on  the  per¬ 
formance  of  the  helicopter.  It  not  only  alters  the  effective 
pitch  angle  of  the  blades  and  thus  the  airloads,  but  also 

Copyright  ©  1989  by  the  American  Institute  of  Aeronautics  and 
Astronautics,  Inc.  No  copyright  is  asserted  in  the  United  States  under  Title 
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produces  highly  nonlinear  interactions  of  the  vortex  with 
the  rotor  flowfield.  It  is  possible  that  such  interactions 
might  produce  vortex  induced  boundary  layer  separation 
which  results  in  loss  of  lift  and  increased  drag.  An  accurate 
simulation  of  the  rotor  flowfield,  therefore,  must  consider 
the  induced  effects  of  the  vortex  wake  including  effects  of 
blade-vortex  interactions  (BVI). 

Numerical  simulations  of  vortex  wakes  is  being  at¬ 
tempted  only  recently  after  bigger  and  faster  computers 
have  become  available  (see  for  example  Ref.  1).  These 
investigations  have  some  limited  success  to  date.  On  the 
other  hand,  much  of  the  progress  in  modeling  the  blade- 
vortex  interaction  has  been  hampered  by  the  lack  of  de¬ 
velopment  of  theoretical  and/or  numerical  techniques  to 
preserve  the  structure  of  the  concentrated  vortices  in  the 
flowfield  without  significant  diffusion.  The  study  of  blade- 
vortex  interaction  has  been  the  subject  of  numerous  re¬ 
cent  research  papers.4-11  These  studies  have  considered 
different  methods  of  preserving  the  concentrated  vortex 
in  the  flowfield.  Among  these,  the  vortex-fitting  technique 
(also  called  the  split  potential  formulation  in  potential  flow 
methods)  has  demonstrated  to  be  a  very  effective  method 
in  preserving  the  vortex  even  when  the  grid  density  is 
sparse.  This  method  has  also  proven  very  economical  7,4 
compared  to  a  more  exact  formulation  used  in  Ref.  5. 

The  purpose  of  this  study  is  to  devise  numerical 
methods  for  the  solution  of  the  Euler  and  Navier-Stokes 
equations  to  accurately  calculate  the  flowfield  of  a  heli¬ 
copter  rotor,  including  the  effects  of  blade-vortex  interac¬ 
tions.  In  particular,  this  paper  concentrates  on  demon¬ 
strating  the  ability  to  calculate  such  an  interaction  flow- 
field  for  parallel  and  oblique  blade- vortex  interactions, 
which  is  three-dimensional  and  unsteady,  on  a  model  heli¬ 
copter  rotor  tested  in  a  wind  tunnel  by  Caradonna  et  al.2'3 
Because  these  experiments  were  specifically  conducted  to 
provide  test  data  to  evaluate  different  numerical  methods, 
the  experimental  apparatus  was  kept  simple  and  consisted 
of  rigid  two-bladed  rotor.  The  blades  are  made  up  of 
NACA  0012  airfoil  and  have  rectangular  planform  with 
no  twist  or  taper;  the  interacting  vortex  was  generated  up¬ 
stream  of  the  model  rotor  by  means  of  a  lifting  rectangular 
wing  mounted  on  the  roof  of  the  tunnel  as  shown  in 
Fig.  2. 

Governing  Equations  and  Numerical  Scheme 

The  governing  partial  differential  equations  are  the 
unsteady,  Euler/thin  layer  Navier-Stokes  equations.  In 
this  study,  the  thin  layer  Navier-Stokes  code  recently  writ¬ 
ten  specifically  for  helicopter  rotor  applications12  will  pro¬ 
vide  the  starting  point  for  further  modifications  to  cal¬ 
culate  the  blade-vortex  interaction  flowfield.  In  order  to 


assess  the  importance  of  viscous  effects  in  these  interact¬ 
ing  flowfields,  both  the  Euler  and  thin  layer  Navier-Stokes 
equations  are  solved.  The  Euler  calculations  are  performed 
by  turning  off  the  viscous  terms  in  the  Navier-Stokes  equa¬ 
tions  set  and  appropriately  modifying  the  surface  bound¬ 
ary  condition.  The  mathematical  formulation  presented  in 
the  following  paragraphs  is  for  the  solution  of  thin  layer 
Navier-stokes  equations  with  the  added  feature  of  intro¬ 
ducing  a  prescribed  vortex. 

For  generality,  the  equations  are  transformed  from 
the  inertial  Cartesian  reference  frame  (*,  y,  z,  t )  to  the  ar¬ 
bitrary  curvilinear  space  ({,  17,  £,  r)  moving  with  the  blade, 
while  retaining  strong  conservation  law-form  to  capture 
shock  waves.  The  transformed  equations,  following  Ref. 
13,  can  be  written  as 

dr (Q  -  Qo)  +  di(F  -  Fo)  +  0„(G  -  Go)  +  d((H  -  Ho) 

=  J ?e-10<S 

(1) 

where 


■  p  ■ 

r  pu 

pu 

,  F  =  J-1 

puV  +  £*p 

pv 

pvU  +  tyP 

pw 

pwU  +  £xp 

.  e  . 

,U(e  +  p)  —  {tp. 

pV  1 

PW 

puV  +  i)xp 

puW  + 

pvV  +  rjyp 

B  -  J- 1 

pvW  +  £yp 

pwV  + 1 jzp 

pwW  +  £,p 

.V(e  +  p)-Tjtp. 

W(e  +  p)  —  £tpJ 

Here  Q  is  the  unknown  flowfield  vector  and  Q„  is  the  so¬ 
lution  of  the  Euler  equations  for  a  presdbed  line  vortex 
aligned  with  the  uniform  free  stream  of  Mach  number, 
Moo j  in  the  y-direction  (see  Fig.  2).  The  terms  F,  G ,  H 
and  F0,  Go,  H0  represent  the  appropriate  flux  vectors  for 
the  two  flows,  respectively.  The  hat,  represented  by  (  ),  in 
Eq.  (1)  denotes  quantities  scaled  by  the  Jacobian  J,  e.g., 
Q  =  J~lQ  etc.  The  contravariant  velocity  components 
U,  V  and  W  are  defined  as 

U  =  6  +  (,u  +  £yv  +  (xw 

V  =TJt+T)oU  +  T)yV  +  T),W 
W  =  Ct  +  C*U  +  CyV  +  (,U> 

In  the  present  formulation,  (  lies  in  the  spanwise  direction, 
1}  is  in  the  chordwise  or  wraparound  direction,  and  £  is 
normal  to  the  blade  surface. 
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The  viscous  flux  vector  5  is  used  in  the  limit  of  thin- 
layer  approximation  and  is  given  by 


S  =  J~l 


0 

+  (p/Z)m2C,x 
/imiVf  +  (/i/3 

+  (p/3)m2(z 

.pm1mi  +  (p./3)m2(£xu  +  <„t>  +  <jU>)- 


where 


the  inertial  frame  of  reference  with  the  boundary  condi¬ 
tions  applied  on  the  rotating  blade.  The  inertial  coordi¬ 
nates  X  =  (x,y,  z,£)  are  related  to  the  blade-fixed  coordi¬ 
nates  Xb  =  (x,y,  z)  through  the  relation  given  by 


X(x,y,z,t)  =  R(t)Xb(x,y,i) 
t  =  i 


(3) 


where  R(t)  is  the  rotational  matrix14  given  by 


mi  =  <1  +  Cy  +  <1 

m2  =  +  (zw( 

m3  =  +  V2  +  w*)(  +  /cPr“1(7  -  l)_1(a2){ 

A 

In  this  equation,  a  is  the  speed  of  sound,  it  is  the  coefficient 
of  thermal  conductivity  and  fi  is  the  viscosity  computed  as 
the  sum  of  laminar  and  turbulent  contributions  for  a  tur¬ 
bulent  boundary  layer.  Sutherland’s  equation  is  used  for 
laminar  viscosity  and  the  two-layer  algebraic  eddy  viscos¬ 
ity  model  of  Baldwin  and  Lomax19  is  used  to  calculate 
the  turbulent  eddy  viscosity,  after  minor  modifications  to 
account  for  the  moving  blade  surface.  Also,  J  is  the  trans¬ 
formation  Jacobian,  whose  inverse  is  written  as 


J  1  =  Hy*\z<.  +  xCJ(zn  +  *>»»<*<  -  xmzn 

~Xvy(Z<  -  X<Vr,Z( 


The  velocity  components  u,  v,  w  and  the  pressure,  p, 
are  related  to  the  total  energy  per  unit  volume,  e,  through 
the  equation  of  state  for  a  perfect  gas  by 


p  =  (7-l)(e-|(«J-MJ +«,’))  (2) 


The  primitive  variables  of  Eq.  (1),  viz.,  the  density 
p,  the  mass  fluxes  pu,pv,pw  and  the  energy  per  unit  vol¬ 
ume  e,  are  normalized  by  the  free  stream  reference  quan¬ 
tities  and  the  pressure  p  by  The  reference  length 

and  velocity  scales  are  the  chord  of  the  blade,  C,  and  free 
stream  speed  of  sound,  ax,  respectively.  The  details  of 
the  metrics  of  transformation  ((« (qt,q*,qy,q*) 
and  (£t,C*,Cy,C*)  be  found  in  Ref.  13. 

In  the  above  equations  u,v,  and  w  are  the  Cartesian 
components  of  velocity  in  the  inertial  coordinate  system 
(x,y,  z,t).  In  the  present  formulation  Eq.  (1)  is  solved  in 


R(t)  = 


cos  fit 
sin  fit 
0 


—  sin  fit  0 
cos  fit  0 

0  1 


(4) 


Here  fl  is  the  angular  velocity  of  the  rotor  and  fit  repre¬ 
sents  the  azimuth  sweep  of  the  rotor  blade. 

The  equations  set,  Eq.  (1),  is  solved  using  a  two- 
factor,  implicit,  finite-difference  numerical  scheme,  devel¬ 
oped  by  Ying  et  al.,ls  that  uses  spatial  central-differencing 
in  the  q  and  (  directions  and  upwind-differencing  in  the  ( 
direction  and  is  given  by 


[/  +  +  hScdn  -  hRe~16cJ-iMnJ  -  £>i|(] 

x  [/  +  h6'(A~)n  +  h6vBn  -  AU]  (A Qn  -  AQB)  = 

-  Af{6|[(P+)"  -  F+]  +  6'[(F-)n  -  F~] 

+  Sn(Gn  -  Go)  +  -  Ho)  -  Re-'S^S?)} 

-  (z>«  I, +  z>,  !<)($"-&) 

(5) 

where  h  =  At  for  first-order  time  accuracy.  Here  6  is 
typically  a  three-point,  second-order-accurate,  central  dif¬ 
ference  operator,  6  is  a  midpoint  operator  used  with  the 
viscous  terms,  and  the  operators  and  6^  are  back¬ 
ward  and  forward  three-point  difference  operators.  The 
time  index  is  denoted  by  n  such  that  t  =  (nAt),  and 
A Qn  —  Qn+1  —  Qn.  The  flux  vector  F  has  been  split 
into  F+  and  F~,  according  to  its  eigenvalues,1®  and  the 
matrices  A ±,  B,  C,  and  M  result  from  the  local  lineariza¬ 
tion  of  fluxes  about  the  previous  time  level  (see  Ref.  13). 
The  finite- difference  scheme  described  in  Eq.  (5)  uses  flux¬ 
splitting  in  (  direction  and  central  differencing  in  the  q  and 
£  directions.  As  a  consequence,  numerical  disspation  terms 
denoted  by  and  Dt  are  employed  in  the  q  and  C,  direc¬ 
tions,  and  are  given  as  combinations  of  second-  and  fourth- 
differences.  In  the  vicinity  of  shock  waves,  the  fourth- 
difference  terms  can  cause  oscillations,  so  it  is  desirable 
to  drop  these  terms  and  rely  only  on  the  second-difference 
terms. 
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The  factored  operators  are  solved  by  sweeping  in  the 
£  direction  and  inverting  tridiagonal  matrices  with  5x5 
blocks  for  the  other  two  directions.  Currently,  a  significant 
part  of  the  computational  time  is  taken  to  form  the  plus 
and  minus  Jacobian  matrices  for  the  flux  vector  F  with  this 
numerical  scheme.  However,  this  effort  has  been  reduced 
by  computing  A+  and  A~  at  every  other  point  (in  both  rj 
and  (  directions)  and  averaging  to  obtain  the  matrices  at 
the  intermediate  points.  The  numerical  code  is  vectorized 
for  the  Cray-2  supercomputer. 

In  writing  Eq.  (5)  it  is  assumed  that  Aj  ~  A±, 
B0  ~  B,  and  C0  ~  C  where  the  Jacobian  matrices 
A0±,  B0  and  C0  correspond  to  the  prescribed-vortex  flow- 
field.  In  the  absence  of  vortex  interaction,  the  prescribed- 
disturbance  flowfield  reduces  to  a  free  stream. 

In  using  this  numerical  algorithm  for  the  physical 
problem  presented  here,  a  note  of  caution  is  in  order.  As 
noted  before,  the  numerical  scheme  is  second-order  accu¬ 
rate  in  space.  Although  it  is  desirable  to  have  higher-order 
accuracy  in  space,  it  turns  out  that  the  present  algorithm 
is  numerically  more  stable  if  this  accuracy  is  decreased 
slightly  from  the  second-order.  Also,  in  evaluating  the 
right  hand  side  of  Eq.  (5),  it  is  desirable  to  calculate  the 
derivatives  of  individual  fluxes  first  and  then  taking  a  dif¬ 
ference.  This  is  particularly  true  for  the  fluxes  in  the  ( 
direction. 

A  body-conforming  finite-difference  grid  hits  been 
used  for  the  rectangular  blade  having  an  aspect  ratio  of  7 
and  consists  of  a  warped  spherical  0-0  grid  topology.  The 
flowfield  grid  is  generated  using  the  three-dimensional  hy¬ 
perbolic  grid  solver  of  Steger  and  Chaussee17  with  proper 
clustering  in  the  leading  and  trailing  edge  regions  and  in 
the  tip  region.  The  grid  is  nearly  orthogonal  at  the  surface 
and  the  spacing  in  the  normal  direction  at  the  surface  is 
chosen  to  be  0.0002C  for  the  viscous  grid  and  0.02C  for  the 
Euler  grid.  The  viscous  computations  are  done  on  a  grid 
consisting  of  21  points  along  the  span,  101  in  the  periodic 
direction  around  the  blade  and  31  in  the  normal  direction. 
The  Euler  grid  is  more  coarser  with  only  21x101x15  grid 
points.  The  outer  grid  boundary  was  chosen  to  be  at  8 
chords  in  all  directions  for  both  viscous  and  inviscid  grids. 

The  boundary  conditions,  both  sin  face  and  farfield, 
are  applied  explicitly.  At  the  rotating  blade  surface  be¬ 
cause  of  the  noslip  condition  for  viscous  flows,  the  con- 
travariant  velocities  U,  V  and  W  are  set  to  zero.  However, 
the  surface  velocity  in  the  inertial  frame  is  nonzero  since  it 
is  equal  to  the  surface  grid  velocity,  determined  by  the  time 
metrics  (t,  rjt  and  as  the  blade  (and  the  grid  attached 
to  it)  moves  in  azimuth.  The  Euler  calculations  use  an  ex¬ 
trapolation  of  contravariant  velocities  to  the  surface.  The 


density  at  the  wall  is  determined  by  zeroth-order  extrapo¬ 
lation.  The  pressure  along  the  body  surface  is  calculated 
from  the  normal  momentum  relation  (see,  for  example 
Ref.  13).  Having  calculated  the  density  and  pressure,  the 
total  energy  is  determined  from  the  equation  of  state. 

At  the  farfield  boundary  the  flow  quantities  are  ei¬ 
ther  fixed  or  extrapolated  from  the  interior  depending  on 
whether  the  flow  is  subsonic  or  supersonic  and  if  it  is  of 
inflow-  or  outflow-type  at  the  boundary.  The  characteristic 
velocities  of  the  Euler  equations  determine  the  number  of 
flow  properties  to  be  specified  to  control  the  reflections  of 
waves  from  the  boundaries.  For  the  subsonic-inflow  bound¬ 
ary,  four  quantities  must  be  specified  and  one  quantity  is 
determined.  The  four  specified  here  are  a  Riemann  invari¬ 
ant,  the  entropy,  and  two  tangential  velocities;  the  quan¬ 
tity  that  is  calculated  is  also  a  Riemann  invariant.  For  the 
supersonic-inflow,  all  flow  quantities  are  specified.  At  the 
subsonic-outflow  boundaries,  only  one  quantity  is  specified. 
For  the  supersonic-outflow  condition  all  flow  quantities  are 
extrapolated  from  the  interior.  The  plane  containing  the 
blade  root  is  chosen  very  close  to  the  rotation  axis  of  the 
blade  (at  R  =  1.0f7)  and  is  also  treated  as  a  farfield  bound¬ 
ary. 

The  interacting  line  vortex  generated  upstream  of  the 
rotor  by  a  rectangular  wing  is  fixed  in  the  inertial  space 
along  a  line  of  constant  x.  It  is  assumed  to  have  an  analyt¬ 
ical  representation  for  the  cylindrical  velocity  distribution 
given  by 


The  constant  Ck  in  the  above  expression  has  been  deter¬ 
mined  to  be  equal  to  0.8  by  matching  the  peak  tangen¬ 
tial  velocity  with  the  experimentally  measured  value  at 
the  measured  radial  distance  reported  in  Ref.  2  by  Srini- 
vasan  et  al.®  The  constant  a0  is  approximately  equal  to  the 
radius  of  the  viscous  core  of  the  vortex  (equal  to  0.167C 
in  the  present  case).  The  induced  velocity  due  to  this  line 
vortex  is  calculated  using  the  Biot-Savart  Law  and  the 
pressure  field  is  calculated  by  solving  the  radial  momen¬ 
tum  equation;  the  vortex  is  initialized  well  upstream  of 
the  interaction  region  as  described  in  Ref.  7. 

Results  and  Discussion 

All  the  calculations  performed  in  this  study  are  done 
in  a  time-accurate  fashion.  The  test  cases  calculated  corre¬ 
spond  to  the  two  sets  of  experiments  of  Caradonna  et  al.2,5 
for  parallel  and  oblique  blade-vortex  interactions.  Since 
one  of  the  purposes  of  the  experiments  was  to  collect  data 
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to  validate  numerical  methods,  the  experimental  appara¬ 
tus  was  kept  simple  to  ease  the  representation  by  numeri¬ 
cal  methods.  The  rotor  geometry  consisted  of  a  two-bladed 
rigid  rotor  of  approximately  14-chord  diameter.  The  blades 
have  a  rectangular  planform  and  consist  of  an  NACA  0012 
airfoil  with  no  twist  or  taper.  The  rotor  blades  are  essen¬ 
tially  nonlifting  in  the  absence  of  the  vortex  interaction. 
The  interacting  vortex  was  generated  upstream  of  this  ro¬ 
tor  by  a  lifting  rectangular  wing  with  an  NACA  0015  air¬ 
foil.  Parallel  and  oblique  BVI  are  generated  by  positioning 
the  line  vortex  along  the  y-axis  (x0  =  0.0)  and  ahead  of 
the  y-axis  (x0  <  0.0)  so  that  the  rotor  blade  encounters 
the  vortex  on  the  advancing  side  in  the  first-  and  second- 
quadrants  of  the  azimuthal  travel  as  shown  in  Fig.  2. 

As  mentioned  before,  the  Navier-Stokes  hover  code 
described  in  Eef.  12  provides  the  starting  point  for  the 
present  investigation.  This  code  has  been  validated  against 
the  experimental  data  of  Caradonna  and  Tung18  for  the 
nonlifting  and  lifting  hover  conditions  at  both  subcritical 
and  supercritical  tip  speeds.  Although  the  main  focus  of 
the  present  investigation  is  to  calculate  both  parallel  and 
oblique  blade-vortex  interactions,  it  is  necessary  to  first 
calculate  the  base-line  rotor-alone  flowfield  solutions  (in 
the  absence  of  vortex  interaction)  at  the  same  free  stream 
conditions.  Two  sets  of  results  for  parallel  blade-vortex 
interactions,  corresponding  to  Mtip  =  0.6  and  0.8  and 
ft  =  0.2  and  one  set  of  results  for  oblique  blade-vortex 
interaction  corresponding  to  Mtip  =  0.763  and  /z  =  0.197, 
will  be  presented  in  the  following  sections  and  compared 
with  experimental  data.  At  the  Reynolds  numbers  cor¬ 
responding  to  the  blade  tip  speeds  in  these  experiments, 
the  boundary  layer  on  the  entire  rotor  blade  could  be  as¬ 
sumed  to  be  turbulent  and,  as  mentioned  before,  an  al¬ 
gebraic  eddy  viscosity  model  of  Baldwin  and  Lomax1*  is 
used  for  calculating  the  turbulent  viscosity. 

Rotor- Alone  Case 

The  objective  of  calculating  the  rotor-alone  flowfield 
solutions  serves  two  purposes.  First,  it  enables  an  under¬ 
standing  of  the  importance  of  unsteady  time-lag  effects 
in  shock  wave  growth  and  decay  as  well  as  the  three- 
dimensionality  of  the  flowfield  of  the  advancing  rotor,  par¬ 
ticularly  when  the  flow  is  supercritical.  Second,  it  provides 
the  baseline  solution  to  eventually  calculate  the  vortex  in¬ 
teraction  flowfield. 

Figure  3  shows  the  instantaneous  surface  pressure 
distributions  for  various  azimuthal  positions  of  the  rotor 
for  the  case  of  Mtip  =  0.6  and  fi  =  0.2.  Examination  of 
the  results  suggests  that  at  this  reference  station  of 
tb  =  0.893  and  for  the  subcritical  flow  condition,  the  flow 
behaves  as  if  it  is  quasi-steady  and  quasi- two-dimensional. 


In  spite  of  the  gradual  increase  of  the  blade-element  rel¬ 
ative  speed  in  the  first  quadrant  and  gradual  decrease  of 
the  same  in  the  second  quadrant,  the  flowfield  appears  to 
remain  nearly  the  same  at  all  azimuth  locations  (as  seen 
in  Fig.  3b)  indicating  that  the  unsteady  time-lag  effects 
are  virtually  absent  for  this  flow  condition  and  the  flow 
behaves  as  if  it  is  quasi-steady.  In  fact,  also  shown  in 
Fig.  3b  is  the  quasi-steady  solution  for  one  azimuth  loca¬ 
tion  of  ip  =  90.°  which  is  in  perfect  agreement  with  the 
unsteady  solutions  at  different  azimuthal  positions.  Also, 
the  comparison  of  three-dimensional  results  of  Fig.  3b  with 
the  two-dimensional  results*  of  Fig.  3a  shows  that  these 
results  are  nearly  identical  suggesting  that  the  flow  also 
behaves  as  if  it  is  quasi -two- dimensional. 

At  the  supercritical  tip  flow  conditions,  correspond¬ 
ing  to  Mtip  —  0.8  and  fj.  =  0.2,  the  basic  rotor  flow- 
field  is  dominated  by  the  presence  of  a  strong  shock  wave 
on  the  advancing  blade  over  large  parts  of  the  first-  and 
second-quadrants.  For  this  advancing  rotor,  the  Euler  and 
Navier-Stokes  results  of  instantaneous  surface  pressures 
are  shown  in  Fig.  4  at  the  radial  station  rg  =  0.893  for 
different  azimuthal  positions  of  the  blade.  The  Euler  re¬ 
sults  show  good  agreement  with  the  experiments  and  the 
Navier-Stokes  results  have  only  a  qualitative  agreement 
because  of  coarse  grid.  It  appears  that  the  shock  wake 
position  is  under  predicted  for  the  Navier-Stokes  method 
and  slightly  over  predicted  for  the  Euler  method.  For 
this  coarse  grid  solution,  the  Navier-Stokes  method  pre¬ 
dicts  the  decay  of  the  shock  much  faster  in  the  second- 
quadrant.  On  the  other  hand,  the  Euler  method  predicts 
that  the  shock  wave  persists  at  the  ip  =  180°  azimuth. 
Previous  two-dimensional  Navier-Stokes  calculations  of  the 
same  flow*  had  shown  strong  three-dimensionality  and  un¬ 
steady  time-lags  in  shock  wave  growth  and  decay.  The 
two-dimensional  assumption  for  this  flow  essentially  over¬ 
predicted  the  shock  wave  position  and  strength,  unlike 
the  subcritical  flow  condition.  In  contrast,  the  three- 
dimensional  Euler  results  seem  to  follow  the  experimental 
observation  correctly. 

Thus  at  this  flow  condition,  the  flowfield  is  highly 
three-dimensional  exhibiting  strong  unsteady  time-lag  be¬ 
haviour  in  the  shock  wave  formation  and  eventual  demise. 
The  unsteady  time-lag  feature  of  the  flowfield  is  demon¬ 
strated  in  the  surface  pressure  distributions  of  Fig.  4 
shown  for  ip  =  60°  and  'P  =  120°  azimuthal  positions, 
which  shows  that  the  two  results  are  very  different  from 
each  other. 

In  principle,  the  shock  wave  should  attain  its  maxi¬ 
mum  strength  when  the  relative  flow  speed  reaches  a  maxi¬ 
mum  value  for  ip  =  90°  blade  azimuthal  position.  However, 
the  shock  wave  keeps  growing  in  strength  even  after  the 
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blade  has  passed  the  ip  =  90°  azimuthal  position  and  a 
weak  shock  wave  seems  to  persist  even  at  the  ip  =  180° 
azimuth.  The  strengthening  of  the  shock  wave  in  the  first- 
quadrant  and  beyond,  and  the  slow  demise  of  it  in  the 
second-quadrant,  suggests  the  existence  of  a  strong  un¬ 
steady  time-lag  effect  at  this  flow  condition.  Such  a  be¬ 
havior  was  totally  absent  at  the  subcritical  flow  condition. 
The  presence  of  the  shock  wave  seems  to  introduce  not 
only  the  time-lag  in  the  adjustment  of  the  flow  as  the  blade 
sweeps  in  azimuth,  but  also  makes  the  flow  highly  three- 
dimensional. 

In  contrast  to  the  flow  on  the  advancing  side  of  the 
rotor  blade,  the  flowfield  of  the  retreating  side  appears 
really  benign  for  this  nonlifting  rotor.  The  lingering  effects 
of  the  shock  wave,  persisting  at  ip  =  180,°  soon  die  out  as 
the  blade  sweeps  into  the  third-quadrant.  Shown  in  Fig.  5 
are  the  surface  pressure  distributions  for  several  azimuthal 
locations  of  the  blade  in  third-  and  fourth-quadrants  on  the 
retreating  cycle.  Since  the  flowfield  is  basically  subcritical 
on  the  retreating  side,  the  plots  of  Fig.  5  appear  to  have 
the  same  behavior  as  those  of  Fig.  3. 

Figure  6  shows  the  surface  pressure  distributions  on 
the  advancing  rotor  at  tb  =  0.946  for  the  flow  conditions 
of  Mtip  =  0.763  and  fi  =  0.197.  Since  this  flow  condition  is 
similar  to  that  of  the  case  of  Fig.  4,  it  has  the  same  quali¬ 
tative  behaviour  of  the  flow  on  the  advancing  side.  As  seen 
in  these  plots,  the  Euler  calculations  are  in  good  agreement 
with  the  experimental  data.3  As  before,  the  Navier-Stokes 
results  have  only  a  qualitative  agreement  for  the  reasons 
mentioned  before.  It  is  important  to  mention  that  the 
rotor-alone  case  involves  a  nonlifting  symmetrical  blade, 
for  which  the  pressure  distributions  should  be  identical  on 
the  upper  and  lower  surfaces.  However,  the  data  in  Fig. 
6,  from  Ref.  3,  seem  to  have  more  scatter  than  the  corre¬ 
sponding  experimental  results  from  Ref.  2. 

Quasi-steady  vs  Unsteady  Results 

Although  the  flowfield  of  a  lifting  helicopter  rotor  is, 
in  general,  unsteady  and  three-dimensional,  a  close  inspec¬ 
tion  of  the  subsonic  surface  pressure  results  presented  in 
Fig.  3b  suggests  that  the  unsteady  and  quasi-steady  re¬ 
sults  are  nearly  the  same  for  this  nonlifting  flow  condition. 
In  fact,  a  quasi-steady  flow  calculated  by  freezing  the  rotor 
blade  at  certain  azimuthal  locations  (say,  ip  =  0.0°  )  pro¬ 
duces  nearly  the  same  flowfield  as  an  unsteady  flowfield  ir¬ 
respective  of  the  flow  conditions.  In  the  present  approach, 
the  unsteady  calculations  are  started  from  a  baseline  quasi¬ 
steady  solution.  At  the  subcritical  flow  conditions,  similar 
to  that  presented  in  Fig.  3,  an  unsteady  calculation  can  be 
started  from  any  azimuthal  location  of  the  blade  because  of 
the  quasi-steady  and  the  unsteady  solutions  being  nearly 


identical  at  all  azimuthal  locations.  However,  this  is  not 
true  for  supercritical  flow  conditions  where  shock  waves 
are  present  in  the  flowfield.  For  example,  Fig.  7  shows  the 
quasi-steady  and  unsteady  surface  pressure  distributions 
for  the  case  of  Mup  =  0.8  and  p  =  0.2  at  azimuthal  loca¬ 
tions  of  ip  =  60°  and  120°.  For  this  flow  condition,  these 
results  indicate  that  the  quasi-steady  assumption  over  pre¬ 
dicts  the  shock  strength  at  ip  =  60°  and  under  predicts  at 
ip  =  120°.  Nevertheless,  at  ip  =  0.0°  position,  the  quasi¬ 
steady  and  unsteady  results  are  nearly  identical  even  for 
this  flow  condition.  So,  for  flow  conditions  that  are  severe, 
it  is  still  reasonable  to  start  the  unsteady  marching  from 
the  quasi-steady  solution  calculated  at  ip  —  0.0°  azimuthal 
location. 

Parallel  Blade- Vortex  Interaction 

During  the  unsteady  three-dimensional  close  en¬ 
counter  of  a  curved  tip  vortex  with  a  rotating  blade,  the 
helicopter  rotor  undergoes  a  variety  of  blade- vortex  inter¬ 
actions  depending  on  the  interaction  angle  between  the 
leading  edge  of  the  rotating  blade  and  the  curved  line  vor¬ 
tex.  These  interactions  are  generally  unsteady  find  three- 
dimensional.  One  limiting  case  of  such  an  encounter,  for 
zero  interaction  angle,  is  termed  parallel  interaction  (see, 
for  example,  Ref.  7).  In  the  experimental  configuration 
considered  here,  this  interaction  occurs  around  ip  =  180° 
azimuthal  position.  For  an  observer  riding  with  the  blade 
(at  a  given  reference  station  along  the  span),  it  appears  as 
though  the  observer  is  passing  a  fixed  vortex  in  the  flow. 
For  of  this  reason,  this  interaction  is  sometimes  approxi¬ 
mated  as  two-dimensional  and  unsteady. 

To  calculate  accurately  the  blade-vortex  interaction 
flowfield,  it  is  necessary  to  preserve  the  vortex  struc¬ 
ture  without  numerical  diffusion.  As  mentioned  before, 
one  method  that  has  been  demonstrated  to  work  effec¬ 
tively  and  economically  in  achieving  this  is  the  prescribed- 
disturbance  scheme.7  The  effectiveness  of  this  scheme  is 
illustrated  in  Fig.  8  (reproduced  from  Ref.  20)  showing 
the  variation  of  lift  coefficient  as  a  function  of  vortex  lo¬ 
cation  during  an  airfoil- vortex  interaction.  Also  shown  in 
this  figure  is  the  Euler  lift  distribution  calculated  using 
a  conventional  (vortex  capturing)  technique.31  Although 
the  two  Euler  solutions  are  computed  on  the  same  grid, 
the  numerical  dissipation  associated  with  the  finite  grid¬ 
spacing  progressively  weakens  the  gradients  and  reduces 
the  effective  vortex  strength  in  the  conventional  method. 
Also,  this  numerical  error  is  grid-dependent  (the  finer  the 
grid  the  lesser  the  error);  however,  it  is  completely  absent 
in  the  prescribed- vortex  solutions,  which  are  essentially  in¬ 
dependent  of  the  grid.  Hence,  the  prescribed-disturbance 
method  has  been  used  here  for  preserving  the  vortex 
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structure. 

Subcritical  Interaction 

The  results  of  a  subcritical  parallel  BVI  are  discussed 
here.  This  case  corresponds  to  the  flow  conditions  of 
Mup  —  0.6,  n  =  0.2,  and  C* T  =  0.133  at  a  blade  refer¬ 
ence  station  (marked  AA  in  Fig.2)  rg  =  0.893.  The  inter¬ 
acting  vortex  is  located  at  x0  =  0.0  and  z0  =  -0.4  along 
the  y-axis.  To  calculate  the  BVI  flowfield,  the  interacting 
vortex  is  initialized  with  the  baseline  rotor  solution  at  the 
azimuthal  position  of  ip  =  120.°  The  unsteady  flow  field  is 
monitored  as  the  blade  continues  to  advance  in  azimuth, 
and  the  peak  interaction  effects  are  shown  in  Fig.  9  in  the 
form  of  instantaneous  surface  pressure  distributions  mea¬ 
sured  around  ip  =  180°  azimuthal  position. 

The  present  Euler  results  are  in  very  good  agree¬ 
ment  with  both  the  experimental  data3  and  a  previous 
two-dimensional  Navier-Stokes  calculation*.  As  seen  from 
Fig.  9a  for  ip  =  178.15,°  the  lift  on  the  blade,  which  is  ini¬ 
tially  zero  in  the  absence  of  vortex  interaction,  is  negative 
(because  of  the  sense  of  rotation  of  vortex)  and  gradually 
becomes  positive  as  the  blade  passes  the  line  vortex.  This 
crossover  of  lift,  from  negative  to  positive,  seems  to  occur 
when  the  vortex  is  approximately  aligned  with  the  quarter- 
chord  line  of  the  blade  (corresponding  io  x„  =  0.0).  As 
mentioned  before,  the  two-dimensional  approximation  of 
this  three-dimensional,  unsteady  interaction  is,  in  fact,  a 
very  good  assumption  for  this  subcritical  flow  condition. 
The  close  agreement  of  the  viscous  and  inviscid  results 
suggests  that  the  viscous  effects  Me  unimportant  for  this 
relatively  weak  interaction. 

Supercritical  Interaction 

This  case  corresponds  to  the  flow  condition  of 
Mtip  =  0.8,  n  =  0.2,  and  C*r  =  0.177  at  a  blade  refer¬ 
ence  station  of  rg  =  0.893.  The  interacting  vortex  loca¬ 
tion  is  same  as  in  subcritical  case.  This  flow  condition,  in 
constrast  to  the  subcritical  flow  condition,  exhibits  strong 
unsteady  time-lag  effects  in  shock  wave  growth  and  decay 
for  the  advancing  rotor.  This  feature  seems  to  accentu¬ 
ate  the  three-dimensional  nature  of  the  flow  and  a  two- 
dimensional  approximation  of  this  flow  over  predicts  the 
shock  wave  strength  and  location  as  demonstrated  in 
Ref.  8. 

As  before,  to  calculate  the  BVI  flowfield,  the  inter¬ 
acting  vortex  is  initialized  at  about  ip  =  120°  azimuthal 
position  and  the  time  evolution  of  the  flow  is  monitored  as 
the  blade  advances  in  azimuth.  Figures  10-11  show  the  in¬ 
teraction  flowfield  results  in  terms  of  instantaneous  surface 
pressures  for  several  azimuthal  positions  of  the  blade  for 
the  Euler  and  the  Navier-Stokes  calculations.  Also  shown 


in  this  figures  are  the  experimental  data3.  Comparison  of 
the  results  shows  that  while  the  Euler  results  (Fig.  10) 
have  good  agreement  with  the  experimental  data  for  all 
azimuthal  locations  of  the  blade,  the  Navier-Stokes  results 
(Fig.  11)  show  only  a  qualitative  agreement  with  the  ex¬ 
perimented  data.  This  feature  is  strongly  dependent  on  the 
accuracy  of  the  base-line  solution  at  the  azimuthal  location 
where  the  vortex  is  initialized.  The  vortex  is  intialized  at 
the  ip  —  120°  position  of  the  blade  for  this  case.  Compar¬ 
ison  of  surface  pressures  of  Fig.  4d  at  this  location  for  the 
Navier-Stokes  and  Euler  calculations  shows  that  the  Eu¬ 
ler  results  have  better  agreement  with  experiments  than 
the  Navier-Stokes  results.  Although  the  Euler  grid  is  very 
coarse,  it  seems  to  capture  the  essential  details  of  the  flow 
sharply  including  the  shock  waves.  Since  the  Navier-Stokes 
grid  is  too  coarse,  the  shocks  appear  smeared.  Neverthe¬ 
less,  it  captures  the  essential  qualitative  features  of  the 
flow. 

The  effect  of  vortex  interaction  is  to  induce  time  de¬ 
pendent  aerodynamic  forces  on  the  rotor.  For  example,  as 
seen  from  the  surface  pressure  plots  of  Figs.  10-11,  the 
lift  on  the  blade,  which  is  intially  zero,  becomes  negative 
(because  of  the  sense  of  rotation  of  vortex)  due  to  induced 
downwash  and  gradually  increases  to  a  positive  value  as 
the  blade  passes  the  fixed  vortex.  The  peak  effects  of  the 
interaction  appear  to  occur  when  the  blade  leading  edge  is 
approximately  above  the  line  vortex. 

Oblique  Blade- Vortex  Interaction 

This  case  corresponds  to  a  free  stream  condition  of 
Mtip  =  0.763  and  p  =  0.197.  The  interacting  vortex  lo¬ 
cated  at  x„  =  -2.13  and  z0  =  -0.25  has  a  strength  of  C*T 
=  0.179  for  the  reference  blade  section  at  rg  =  0.946. 
In  contrast  to  the  parallel  interaction,  the  oblique  inter¬ 
action  occurs  over  a  larger  azimuthal  sweep  of  the  rotor 
blade  starting  in  first-quadrant  and  completing  in  second- 
quadrant.  The  peak  interaction  occurs  at  around  ip  =  20° 
in  the  first-quadrant  and  around  ip  =  160°  in  the  second- 
quadrant.  Therefore,  the  vortex  is  initialized  at  ip  =  0.0° 
and  the  evolution  of  the  unsteady  interacting  flow  is  mon¬ 
itored  for  the  advancing  rotor. 

Figures  12-13  show  the  surface  pressure  distributions 
for  several  azimuthal  positions  of  the  blade  from  the  Euler 
and  the  Navier-Stokes  calculations  of  the  interaction  flow- 
field.  As  observed  before,  the  transonic  Euler  calculations 
agree  better  with  the  experimental  data*  than  the  Navier- 
Stokes  calculations  for  the  reasons  cited  earlier.  On  the 
other  hand,  the  interaction  in  the  first-quadrant  occurs  at 
subcritical  flow  condition  and  therefore  the  agreement  of 
both  the  Euler  and  Navier-Stokes  calculations  with  the  ex¬ 
perimental  data  is  good  at  ip  =  50°.  For  the  interaction 
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in  the  second-qudrant,  where  the  flow  is  supercritical,  the 
Euler  results  have  better  agreement  with  the  experiments 
than  the  Navier-Stokes  results.  Compared  to  the  preceding 
transonic  parallel  interaction  results,  the  Navier-Stokes  re¬ 
sults  for  this  interaction  appear  to  agree  better  with  the 
experimental  data.  This  is  in  direct  consequence  of  the 
better  agreement  of  the  rotor-alone  Euler  and  the  Navier- 
Stokes  surface  pressures  at  all  azimuthal  positions,  as  seen 
in  Fig.  12  for  this  flow  condition.  The  maximum  relative 
blade  tip  speed  in  this  case  (for  ip  =  90°  position)  is  0.91 
compared  to  a  value  of  0.96  corresponding  to  the  transonic 
parallel  interaction  case  and  is  thus  a  weaker  interaction. 

From  the  surface  pressures  of  Fig.  13,  it  can  be  seen 
that  the  lift  on  rotor  blade,  which  is  initially  zero,  be¬ 
comes  negative  due  to  the  vortex  influence  (because  of  the 
sense  of  rotation),  increases  to  a  maximum  (negative)  value 
in  second-quadrant,  and  changes  over  to  a  positive  value 
at  about  ip  =  165°  azimuthal  position  of  the  blade.  The 
peak  interaction  appears  to  occur  around  the  ip  =  160° 
azimuthal  position. 

Conclusions 

A  numerical  procedure  is  presented  to  calculate  the 
unsteady,  viscous  and  inviscid,  three-dimensional  flowflelds 
of  a  helicopter  rotor  in  forward  flight  and  parallel  and 
oblique  blade- vortex  interactions  in  subcritical  and  super¬ 
critical  flow  conditions.  Important  flow  features  such  as 
unsteady  time-lag  effects  in  shock  wave  growth  and  demise, 
as  well  as  the  importance  of  three-dimensional  effects  are 
discussed.  While  it  is  possible,  under  certain  flow  con¬ 
ditions,  to  approximate  the  parallel  blade-vortex  interac¬ 
tion  as  two-dimensional  and  unsteady,  the  oblique  blade- 
vortex  interaction,  on  the  other  hand,  is  naturally  three- 
dimensional,  and  unsteady  and  has  to  be  treated  that  way 
without  any  simplifications. 

The  numerical  results  are  compared  with  two-sets  of 
test  data  generated  by  Caradonnaet  aLJ,s  on  a  model  two- 
bladed  rotor  in  a  wind  tunnel.  While  the  Euler  results 
have  good  agreement  with  experiments,  the  Navier-Stokes 
results,  with  a  coarse  grid,  have  only  a  good  qualitative 
agreement.  Nevertheless,  the  study  demonstates  that  the 
numerical  methodology  presented  here  is  capable  of  calcu¬ 
lating  accurately  the  flowfield  of  a  helicopter  rotor  includ¬ 
ing  the  effects  of  vortex  interaction  under  subcritical  and 
supercritical  flow  conditions. 
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Fig.  1  Schematic  of  the  complex  flowfield  of  a  helicopter  rotor  in  forward  flight. 
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Fig.  3  Instantaneous  surface  pressure  distributions  at  the  reference  blade  section  for 
advancing  rotor.  Mtip  =  0.6,  p  =  0.2,  Re  =  2.1x10*,  rB  =  0.893. 
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- navier-stokes  calculations 

- EULER  CALCULATIONS 

•  EXPERIMENT  (CARADONNA  et  al.) 
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Fig.  4  Instantaneous  sui '«.„■’  pressure  distributions  at  the  reference  blade  section  for 
advancing  rotor.  Mtv  —  0.8,  fi  =  0.2,  Re  =  2.89x10®,  rs  =  0.893. 
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Fig.  5  Instantaneous  surface  p 
retreating  blade-Euler  solution. 


-  NAVIER-STOKES  CALCULATIONS 

-  EULER  CALCULATIONS 

•  EXPERIMENT  (CARADONNA  et  at.) 


Fig.  6  Instantaneous  surface  pressure  distributions  at  the  reference  blade  section  for 
advancing  rotor.  M«p  =  0.763,  p  =  0.197,  Re  =  2.75x10*,  j\b  =  0.946. 
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—  EULER  UNSTEADY 
EULER  QUASISTEADY 
•  EXPERIMENT  (CARADONNA.  at  al.) 
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Pig.  7  Comparison  of  pressure  distributions  for  quasi-steady  and  unsteady  Euler 
calculations.  M«p  =  0.8,  (i  =  0.2,  tb  =  0.893. 


- EULER  (PERTURBATION) 

- EULER (NONPERTURBATION) 
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Fig.  8  Effectiveness  of  the  presdbed^vortex  method  for  two-dimensional  BVI: 

NACA  64A006  airfoil.  Mx>  =  0.85,  T  =  0.2,  z0  =  -0.26.  (Reproduced  from  Ref.  20.) 
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—  NAVIER-STOKES  RESULTS  -  PARALLEL  INTERACTION 


Fig.  11  Instantaneous  surface  pressures  during  parallel  blade- vortex  interaction  -  the 
Navier-Stokes  results.  The  solid  line  and  the  solid  circle  refer  to  the  vortex  side  of  the 
blade.  Mti,  =  0.8,  n  =  0.2,  Re  =  2.89x10*,  f  =  0.177,  x0  =  0.0,  z,  =  -0.4,  rB  =  0.893. 
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-  ---  NAV1ER-STOKES  RESULTS  -  OBLIQUE  INTERACTION 

•  o  EXPERIMENTS  (CARADONNA  et  al.) 


Fig.  13  Instantaneous  surface  pressures  during  oblique  blade-vortex  interaction  -  the 
Navier-Stokes  results.  The  solid  line  and  the  solid  circle  refer  to  the  vortex  side  of  the 
blade.  Mti,  =  0.763,  n  =  0.197,  Re  =  2.75x10*,  f  =  0.179,  c.  =  -2.13,  r.  =  -0.25, 
ra  =  0.946. 
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